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A  SYSTOLIC  2-D  CONVOLUTION  CHIP 

H.  T.  Rung  and  S.  W.  Song 
Department  of  Computer  Science  * 

Carncgie-Mellon  University 
Pittsburgh,  Pa.  15213 

Last  Revised  March  1981  , 

Abstract: 

This  paper  describes  a  chip  for  performing  the  2-D  (two-dimensional)  convolution  in  signal  and  image 
processing.  The  chip,  based  on  a  systolic  design,  consists  of  essentially  only  one  type  of  simple  cells,  which  are 
mesh-imerconnectcd  in  a  regular  and  modular  way,  and  achieves  high  performance  through  extensive 
concurrent  and  pipelined  use  of  these  cells.  Denoting  by  u  the  cycle  time  of  the  basic  cell,  the  chip  allows 
convolving  a  kxk  window  with  an  nxn  image  in  0(n2u/k)  time,  using  a  total  of  k3  basic  cells.  The  total 
number  of  cells  is  optimal  in  the  sense  that  the  usual  sequential  algorithm  takes  0(n2k2u)  time.  Furthermore, 
because  of  the  modularity  of  the  design,  the  number  of  cells  used  by  the  chip  can  be  easily  adjusted  to  achieve 
any  desirable  balance  between  I/O  and  computation  speeds. 

Keywords:  * 

Image  and  signal  processing,  parallel  computation,  special-purpose  chip  design,  systolic  arrays,  two- 
dimensional  convolution,  VLSI  designs. 

1 .  Introduction 

With  recent  technological  advances  in  the  VLSI  circuitry,  the  chip  capacity  (or  component  count  on  a  chip) 
is  increasing  at  an  astonishing  rate  [7],  Doth  the  opportunities  and  challenges  regarding  effective  use  of  VLSI 
are  tremendous.  Systolic  design  is  an  architectural  concept  proposed  for  VLSI  (3).  Chip  designs  based  on 
systolic  architectures  tend  to  be  simple,  modular,  and  of  high  performance.  In  [6j  a  general  discussion  on  the 
attractiveness  of  systolic  architectures  is  given.  Systolic  architectures  are  particularly  suited  to  chip 
implementation  of  operations  in  signal  and  image  processing  such  as  filtering,  correlation,  and  discrete 
Fourier  transform  (5].  In  this  paper  we  introduce  a  chip,  based  on  a  novel  systolic  design,  for  performing  the 
2-D  convolution  operator.  Examples  of  application  of  this  operator  include  noise  smoothing,  linear  edge 
enhancement,  edge  crispcning,  etc..  Similar  designs  can  be  used  for  digital  filtering  and  numerical  relaxation. 
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2.  Description  of  the  Systolic  Design 

To  define  the  2-D  convolution  problem,  consider  a  matrix  (or  image)  X  =  {x^}  and  a  kxk  window  with 
weighting  coefficients  as  shown  in  Figure  2- 1.  For  the  purpose  of  illustration,  we  use  k  =  3.  Now  slide  the 
3x3  window  along  the  matrix.  For  each  position  of  the  window,  the  weighted  sum  of  the  entries  (or  pixels)  in 
the  submatrix  covered  by  the  window  is  to  be  computed.  More  precisely,  if  the  entry  at  the  center  of  the 
submatrix  is  x  ,  then  we  wish  to  compute 

2,3=1  2j3=j  «ij  xr+-i-2_  s+j-T 


X 

X 

X 

•  •  .  X 

'll 

12 

13 

14 

In 

X 

X 

X 

.  .  .  X 

21 

22 

23 

24 

2 n 

X 

X 

X 

.  .  .  X 

31 

32 

33 

34 

3n 

X 

X 

X 

.  .  .  X 

nl 

n2 

n3 

n4 

nn 

W11 

U) 

12 

0) 

13 

U21 

22 

W23 

W31 

<a) 

32 

U33 

WEIGHTING 

COEFFICIENTS 


3x3  WINDOW 


Figure  2-1:  nxn  matrix  and  kxk  window  with  k  =  3. 


The  2-D  convolution  chip  is  based  on  a  variant  of  the  systolic  FIR  filtering  array  proposed  in  [4]  and  [5]. 
The  chip  uses  essentially  only  one  type  of  basic  cells  (see  Figure  2-2),  which  are  interconnected  in  a  regular 
and  modular  way  to  form  a  two-dimensional  systolic  array.  The  function  of  a  basic  cell  is  to  update  a  result  Y 
as  indicated  in  the  figure.  Note  that  in  each  cell,  w  is  a  preloadcd  weighting  coefficient  and  each  x  stays  inside 
the  ceil  for  a  cycle. 
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Figure  2-2:  Basic  cell. 
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RESULT 


Figure  2-3:  Systolic  array  consisting  of  three  kernel  cells. 


The  three  kernel  cells  that  form  the  systolic  array  are  shown  in  Figure  2-3.  Each  kernel  cell  is  composed  of 
nine  basic  cells  and  one  row-interface  cell  whose  function  will  become  clear  later.  Five  rows  of  the  matrix  X 
advance  synchronously  from  left  to  right.  The  small  dots  in  the  figure  denote  appropriate  delays  needed  to 
make  each  column  advance  as  a  whole.  (This  does  not  mean  that  these  delays  have  to  be  implemented  on- 
chip:  the  inputs  can  be  spaced  accordingly  when  fed  into  the  device.)  During  the  cycle  subsequent  to  that  an 
entry  enters  the  upper  right  basic  cell  of  a  kernel  cell,  the  weighted  sum  corresponding  to  a  submatrix  with 
that  entry  as  the  upper  left  entry  will  be  output  from  that  kernel  cell.  Therefore  by  sweeping  the  first  five 
rows,  a  systolic  array  consisting  of  three  kernel  cells  can  compute  the  first  three  output  rows  on-the-fly.  In 
general,  by  sweeping  rows  3i+ 1  to  3i+ 5  the  array  computes  rows  3i+ 1  to  3i+ 3.  Thus  an  entry  in  the  nxn 
matrix  is  input  to  the  cell  array  at  most  twice.  In  the  rest  of  the  section,  we  describe  the  underlying  idea  of  the 
kernel  cell  design. 

Suppose  we  want  to  compute  the  following  weighted  sum  in  the  first  row  of  a  kernel  cell: 

Y  =  WU  Xij-2+  WUXij-l  +  W13Xij- 

One  way  to  do  this  is  to  use  a  special  case  of  the  linear  systolic  FIR  filtering  array  discussed  in  [5],  We  use  a 
two-way  flow  basic  cell  (to  distinguish  it  from  the  one-way  flow  basic  cell  defined  earlier)  as  in  Figure  2-4,  and 
the  desired  value  for  Y  can  be  obtained  as  illustrated  in  Figure  2-5  (a)  through  (c).  In  Figure  2-5  (a),  we  have 
denoted  by  Y  the  particular  Ym  to  the  rightmost  basic  cell,  with  its  value  initialized  to  zero. 

An  improved  way  is  to  use  the  one-wav  flow  basic  cell  as  defined  in  Figure  2-2.  Consider  Figure  2-6  (a) 
where  three  basic  cells  in  the  first  row  of  a  kernel  cell  jre  shown.  Consider  the  moment  the  entry  x(.  enters  the 
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Figure  2-4:  A  two-way  flow  basic  cell. 


(a) 


(b) 


(O 


Figure  2-5:  Computing  Y  with  two-way  flow. 

leftmost  basic  cell.  Denote  by  Y  the  particular  Yjn  to  this  ceil,  whose  value  is  initialized  as  zero.  During  the 
cycle  Xg  enters  the  leftmost  cell.  Y  becomes  w,3  x^.  In  the  next  cycle  xjj  3  is  accumulated  to  Y.  as  shown 
in  Figure  2-6  (b).  Finally,  as  in,  Figure  2-6  (c),  x.  2  is  further  accumulated  to  Y.  As  a  result  when  Y 
emerges  from  the  rightmost  cell,  it  has  already  accumulated  the  weighted  sum  of  the  first  row  of  a  submatrix. 
In  the  meanwhile,  the  same  has  been  happening  at  the  two  bottom  rows  of  the  kernel  cell.  Therefore,  by 
summing  the  partial  results  at  the  row-interface  cell,  the  desired  output,  which  is  the  weighted  sum  of  all  the 
nine  entries  in  u  1x3  submatrix,  is  obtained  (see  Figure  2-3). 


The  filtering  array  of  Figure  2-6  differs  from  the  previous  one  in  that  both  the  x-stream  and  Y-stream  now 


Figure  2-6:  Computing  Y  with  one-way  flow. 

move  from  left  to  right  and  the  x-stream  travels  at  half  the  speed  of  the  Y-stream.  This  variation  results  in  the 
active  use  of  every  cell  at  any  given  time.  As  a  result  each  kernel  cell  can  produce  a  pixel  every  cycle  rather 
than  every  other  cycle,  though  each  basic  cell  now  requires  the  additional  storage  to  provide  the  needed  delay 
for  the  x-strcam. 

3.  Balancing  I/O  and  Computation 

Because  of  the  modularity  of  the  basic  design  of  the  systolic  convolution  array,  its  size  can  easily  be 
adjusted  to  match  the  bandwidth  between  the  array  and  the  host  to  which  it  is  attached.  Let  u  denote  the 
cycle  time  of  the  basic  cell.  Suppose  that  the  image  is  nxn,  and  that  in  time  u,  2k- 1  pixel  values  can  be 
transferred  from  the  host  to  the  array,  and  k  output  pixels  can  be  sent  back  to  the  host.  For  processing  the 
whole  image,  the  I/O  time  alone  is  thus  Ofn~ u/k).  To  balance  this,  the  convolution  array  allows  convolving  a 
kxk  window  with  the  whole  image  in  0(n*u/k)  time,  using  k  kxk  kernel  cells.  We  thus  see  how  I/O 
complexity  can  be  used  to  guide  the  design  of  special-purpose  devices.  The  reader  is  referred  to  [2]  for  several 
lower  bound  results  on  the  I/O  complexity  for  a  number  of  problems  including  the  fast  Fourier  transform. 
The  total  number  of  cells,  kj,  used  in  the  proposed  convolution  array  is  optimal  in  the  sense  that  the  usual 
sequential  algorithm  takes  Ofn'k’u)  time. 
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4.  Implementation 

Using  part  of  a  chip  we  have  implemented  a  3x3  kernel  cell,  consisting  of  nine  basic  cells  plus  one  row- 
interface  cell.  The  implementation  uses  a  bit-serial  word-parallel  organization.  Pixels  are  input  as  8-bit 
samples  and  output  with  16  bits.  The  weighting  coefficients  can  be  changed  during  the  loading  phase  prior  to 
the  convolution  computation.  Their  values  are  restricted  to  be  powers  of  two  in  [- 16,  16],  to  reduce  the  area  of 
multiplier  circuits  in  each  basic  cell.  (See  the  remark  in  Section  4.4  for  a  technique  of  using  the  chip  to  handle 
cases  where  weighting  coefficients  arc  not  powers  of  two.) 

We  expect  that  based  on  this  prototype  design  a  full  chip  can  contain  three  3x3  kernel  cells  and  output  a 
pixel  in  less  than  175  ns.  The  anticipated  high  computation  rate  is  mostly  due  to  the  high  degree  of 
concurrency  inherent  in  the  overall  design. 


4.1 .  I/O  Description 

A  total  of  twelve  I/O  pins  are  used.  Eleven  of  these  are  for  input;  only  one  is  for  output.  Weighting 
coefficients  are  loaded  into  the  convolution  chip  prior  to  the  computation  proper  and  therefore  the  inputs  of 
these  coefficients  and  the  pixel  values  share  the  same  pins.  A  cycle  is  composed  of  16  minor  cycles  each  of 
which  includes  phases  <pL  and  <pr  During  each  cycle  one  pixel  value  is  input  and  one  result  is  output.  Since 
an  input  pixel  value  has  only  8  bits  which  are  input  bit  serially,  input  will  occur  every  other  half-cycle.  In  the 
following  we  give  a  list  of  names  of  each  pin  accompanied  by  a  brief  description.  Some  of  the  control  signals 
can  in  fact  be  generated  on-chip. 

Name  Description 

Vdd  Power 

Gnd  Ground 

<pL  Clock  phase  1 

<p2  Clock  phase  2 

xor«1  Pixel  value  or  weighting  coefficient  for  row  1 

xorw2  Pixel  value  or  weighting  coefficient  for  row  2 

xorw3  Pixel  value  or  weighting  coefficient  for  row  3 

Loaded  Control  signal  indicating  weighting  coefficients  are  being  loaded 

LSB  Control  signal  indicating  the  least  significant  pixel  bits  arc  being  input 

Circulate  Control  signal  indicating  the  2nd  1/2  cycle  (absence  of  input) 

Yaux  Optional  input  to  be  accumulated  to  the  computed  result 

(to  allow  repeated  computation  for  non-power-of-2  weights)  < 

Y  Result  of  the  weighted  sum 
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Circulate 


Circulate 


Figure  4-1:  General  layout  of  a  kernel  cell. 

4.3.  Implementation  of  the  Basic  Cell 

The  basic  cell  is  represented  schematically  as  in  Figure  4-2.  Each  pixel  value  remains  in  a  basic  cell  for  two 
cycles  while  each  result  stays  there  for  one  cycle.  Two  8-bit  shift  registers  are  provided  for  the  storage  of  two 
8-bit  pixel  values.  During  the  first  half  cycle,  pixel  bits  entering  a  basic  cell  advance  from  left  to  right  into  the 
first  8-bit  shift  register,  and  at  the  same  time  enter  the  multiplier,  as  indicated  by  the  solid  lines.  Dunng  this 
same  time,  pixel  bits  already  inside  the  first  8-bit  shift  register  are  moved  into  the  second  S-bit  shift  register. 
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Basic  celt 


flow  during  1st  1/2  cycle 
How  during  2nd  1/2  cycle 


flow  during  whole  cycle 


flow  during  the  loading  stage  for  loading 
the  weighting  coefficients 


Figure  4-2:  Implementation  of  the  basic  cell. 

while  pixel  bits  originally  inside  the  second  8-bit  shift  register  are  moved  to  the  next  basic  cell.  During  the 
second  half  cycle  all  pixel  bits  circulate  inside  the  8-bit  shift  registers.  Flow  during  the  second  half  cycle  is 
indicated  by  the  dash  lines. 


Multiplication  in  this  particular  case  is  merely  a  shift  operation.  It  takes  one  minor  cycle  for  an  input  bit  to 
enter  the  shifter  and  the  result  bit  to  exit  the  serial  adder.  A  15-bit  shift  register  therefore  is  prov  ided  to  hold 
the  result  bits  so  that  they  will  arrive  at  the  next  basic  cell  at  the  appropriate  time. 
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4.4.  The  Row-Interface  Celt 

This  cell  computes  the  sum  of  the  three  partial  results  from  the  three  rows  of  basic  cells.  It  also  takes  an 

optional  input  Yaux  which  will  be  accumulated  to  the  computed  result  This  will  allow  the  handling  of  cases 

in  which  weighting  coefficients  are  not  powers  of  two  by  repeating  the  computation  several  times.  In  most 

image  processing  applications  it  is  sufficient  to  have  weighting  coefficients  as  powers  of  two  and.  in  such 

cases.  Y,  will  be  zero, 
aux 

5.  Concluding  Remarks 

We  have  developed  a  chip  design  for  the  2-D  convolution  operator,  which  uses  essentially  only  one  type  of 
basic  cells  interconnected  in  a  regular  and  modular  way.  While  this  experiment  is  about  a  particular  design, 
we  wish  to  make  the  following  general  remark  concerning  the  balance  of  I/O  and  computation.  Since  a 
special-purpose  device  is  typically  attached  to  a  host,  from  which  it  gets  the  data  to  be  processed  and  to  which 
it  outputs  results,  I/O  considerations  play  an  important  role  on  the  overall  performance.  The  ultimate  goal  of 
a  special-purpose  hardware  design  is  (and  should  be  no  more  than)  that  of  achieving  a  computation  rate  that 
balances  the  available  I/O  bandwidth.  Thus  it  is  important  that  the  design  be  modular  so  that  its  size  can 
easily  be  adjusted  to  match  the  bandwidth  between  the  device  and  the  host.  In  the  design  presented  basic 
cells  are  implemented  in  a  bit-serial  manner.  As  a  result,  the  time  fora  basic  cel)  to  process  a  pixel  is  likely  to 
be  longer  than  that  for  the  memory  to  input  and  output  a  pixel.  This  calls  for  parallel  inputs  and  outputs  of 
multiple  pixels  at  each  cycle  of  the  basic  cell.  This  is  the  reason  why  we  made  the  systolic  array  of  this  paper 
input  five  pixels  and  output  three  pixels  every  cycle.  If,  on  the  other  hand,  basic  cells  are  implemented  by  bit- 
parallel  schemes  (which  imply  a  much  shorter  cycle  time  for  the  basic  cell),  then  the  chip  is  likely  able  to  input 
and  output  only  one  pixel  during  each  cycle.  A  one-dimensional,  rather  than  two-dimensional,  systolic  array 
would  be  an  appropriate  design  for  this  case. 
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